o 
o 

(N 



> 

(N 

o 
\o 
o 

o 

:a 

Ph: 



X 



QCD pressure up to four loops at finite temperature and density 
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In recent years the perturbative expansion of the pressure of massless QCD has been driven to order g^lng at 
high temperatures and finite chemical potentials, which has required calculations up to three-loop order in the full 
theory and up to four-loop order in three-dimensional effective theories. In the present paper we briefly review 
the theoretical background behind this work and explain some of the methods used in the computations. 



1. INTRODUCTION 

In order to obtain a quantitative picture of the 
behavior of the QCD pressure in the entire de- 
confined phase — not just in the region /i « 0, 
T ^ Tc that is accessible by lattice simulations 
— one needs to continuously work to drive its 
perturbative expansion to new orders. At high 
temperatures and zero chemical potential it has 
been long known that this expansion is plagued 
by bad convergence properties and a large renor- 
malization scale dependence, but the highest or- 
der [g^\ag) results [1] derived so far have given 
reason for a new hope that the inclusion of the 
next 0{g^) term will improve the situation con- 
siderably [2]. Furthermore, one has been able 
to show that the slow convergence is due to the 
softest bosonic degrees of freedom, which implies 
that purely fermionic observables, such as quark 
number susceptibilities [3,4], behave in a much 
cleaner way and that even a simple reorganiza- 
tion of perturbation theory may have significant 
effects on the properties of the pressure. Hence 
there is an obvious motivation to extend the cal- 
culations to the next order, where the first fully 
non-perturbative contributions from the ultrasoft 
scale g'^T enter the picture. 

Whereas the high-temperature pressure has 
been almost constantly under active research dur- 
ing the last three decades, perturbation theory at 
high densities and small temperatures has drawn 
very little attention since the late 1970's. Apart 
from the emergence of the fundamentally non- 
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perturbative phenomenon of color superconduc- 
tivity in this region of the ^-T plane, this is to a 
large part due to certain computational problems: 
having no effective lower-dimensional theories to 
work with, one is forced to tackle the inevitable 
IR problems by explicit resummations of infinite 
classes of diagrams. This has proven to be a chal- 
lenging numerical task even at order 17'*, and it 
is only very recently that any progress has been 
made at small but non-zero temperature [5] . 

The purpose of the present report is to in- 
troduce to the reader the most important ma- 
chinery of finite-temperature perturbation the- 
ory with special emphasis on the diagrammatic 
tools needed in the determination of the QCD 
pressure at high T. We start by reviewing the 
concept of dimensional reduction and carry on 
to describe the four-loop diagrammatic calcula- 
tions one encounters both in full QCD and in its 
three-dimensional effective theories when trying 
to extend the perturbative series of the pressure 
to order g^ . To conclude, we then briefly outline 
some recent efforts to study the behavior of the 
low-temperature pressure perturbatively. 

1.1. Basics of flnite-T field theory 

The most fundamental quantity that describes 
the equilibrium thermodynamics of a grand 
canonical ensemble is its grand potential 

n{V,T,^li) = -TlnZ, (1) 

where Z denotes the partition function 

Z{V,T,{^i^}) = TTe^p\~p{H-Y,^^^N^)], (2) 



f3 = 1/T and the chemical potentials fii cor- 
respond to conserved charges. In analogy with 
the generating functional of the zero-temperature 
Green's functions, one can easily derive a func- 
tional integral representation for the partition 
function of an interacting quantum field theory 



Z = 



P^exp 



dr 



d^xC 



(3) 



with (/) being periodic on the interval [0, (3] for 
bosons and antiperiodic for fermions. This 
straightforwardly leads to the finite-temperature 
Feynman rules, which differ from the zero-T ones 
only by the fact that the po integrals are replaced 
by discrete sums 



dpo 

27T 



^E 



(4) 



Po 



over the Matsubara frequencies 



{po)bos = 2TmT - ifibos, n eZ (5) 

(po)/er = {2n+l)TTT -ififer, n eZ. (6) 

The partition function is finally available through 
the computation of the one-particle irreducible 
(IPI) vacuum graphs of the theory. 

2. LARGE T/fi: DIMENSIONAL RE- 
DUCTION 

Our eventual goal is to determine the pertur- 
bative expansion of the pressure of hot QCD up 
to order g^ both at zero density and at finite but 
moderate quark chemical potentials /i/ < lOT, 
which we in the following will assume to be the 
case. It has been known for a long time that in 
this region the result is not simply a power series 
in g^ — which one would naively expect — but 
instead contains also odd powers and logarithms 
of the coupling constant 

PQCD = T*[Ao + A2g^ + A3g'' + A'^gHng 
+Aig^ + A5.g' + ^6.9'ln.g + ^6.g' + Oig')] . (7) 

The reason for this is the infrared sensitivity of 
the quantity, which makes its straightforward di- 
agrammatic expansion diverge already at three 
loops. At this order one namely needs to take 



into account the screening of electrostatic gluons 
(for which moebye ~ gT) and at four loops also 
the non-perturbative screening of the magneto- 
static ones {m-fnagn ~ 9^T), which implies the ne- 
cessity to re-organize the perturbative expansion. 
The most natural and straightforward means to 
achieve this is to apply the machinery of dimen- 
sional reduction and effective three-dimensional 
theories, which in a conceptually simple way leads 
to the emergence of the terms non-analytic in g^. 

In its most straightforward formulation dimen- 
sional reduction is simply a statement of the ob- 
vious fact that if the temperature is considerably 
higher than any other energy scale in the sys- 
tem (in our case most importantly Aqcd), the 
non-static modes of the different fields effectively 
decouple because of their large thermal masses. 
This leaves as the dominant degrees of freedom 
of high-temperature QCD the electrostatic (^o) 
and magnetostatic {Ai) gluons, which, as noted 
above, however do not remain massless them- 
selves; in order for the setup to work one needs 
to assume a (parametrically) clear scale hierar- 
chy from the onset: g'^T <^ gT <^ 2ttT. One 
may then successively integrate out the scales 
27tT and gT, leading to a sequence of three theo- 
ries: full QCD describing the hard scales ~ 27rT, 
a three-dimensional (because the temporal direc- 
tion has vanished) Yang-Mills + adjoint Higgs 
theory (EQCD) corresponding to the soft ones 
~ gT, and a 3d pure Yang-Mills theory (MQCD) 
describing the ultrasoft scales ~ g^T. This list is 
exhaustive, since due to confinement there are no 
more softer scales in the system. 

For 11 — the details of the derivation of 
the effective theories and their Lagrangians can 
be found in Refs. [1,6], and the minor modifi- 
cations that finite chemical potentials induce in 
Refs. [4,7]. The implications the procedure has 
for the computation of the full theory pressure are 
the following: it can be written as a sum of three 
terms, pqcd = Pe + Pm + Pg, that correspond 
to the scales 2TrT, gT and g^T, respectively. Be- 
cause all the non-perturbativity of the pressure 
is contained in the last one, the first two can 
be computed as strict perturbative expansions, 
i.e. pure diagrammatic expansions, in the corre- 
sponding theories. Both the IR and UV diver- 
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Figure 1. The ghost-free one-, 
loop graphs contributing to pe ■ 



two- and three- 



gences of these theories can be handled through 
dimensional regularization, because in the even- 
tual sum of the three terms they cancel with each 
other. The function pc, which denotes the pres- 
sure of MQCD and gives its first contributions at 
orders g^hig and g®, has to be computed non- 
perturbatively, either though a combination of 
ordinary lattice simulations and lattice pertur- 
bation theory or through stochastic perturbation 
theory [8]. It is, however, notable that the whole 
contribution of the Lindc sea [9] that for a long 
time was believed to render the entire expansion 
of the pressure useless has now been separated 
into a well-defined and computable (though yet 
unknown) number B2 



PG = ff^[Biln.g + B2]+0(/). 



(8) 



2.1. The hard scales 

As noted above, the function pE is computable 
through the strict perturbation expansion of the 
full theory pressure, which up to three-loop order 
(or, equivalently, g^lng) is given by the graphs 
of Fig. 1. All one needs to do in principle is to 
write down the expressions of the diagrams, per- 
form the color and Lorentz algebras and finally 
evaluate the sum-integrals analytically. At the 
present order this is indeed enough, as — at least 



in the Feynman gauge — the number of indepen- 
dent sum-integrals that remain to be computed 
is low enough to allow for their explicit evalua- 
tion by hand. At three loops these 'master' sum- 
integrals are actually all of the type 



{{p-Qfr 



\2 ' 



(9) 

^PQR P^Q^ {rT (p - RY {Q-RY' 

with P, Q and R being either bosonic or fcrmionic 
and the index m equalling 0, 1 or 2. 

In brief, the general strategy in the analytic 
calculation of sum-integrals of the above type, 
first developed for vanishing chemical potentials 
in Ref. [10] and later generalized to /k ^ in 
Refs. [4,11], amounts to 

• employing Lorentz invariance and elemen- 
tary linear changes of integration variables 
to write the integrands in terms of scalar 
'polarization' functions of the type 



n(p) 



1 



Q^iQ-P? 



(10) 



where no spatial components of the mo- 
menta appear in the numerators, 

• subtracting the leading UV parts of the po- 
larization functions from the sum-integrals 
and treating the corresponding divergent 
terms analytically in momentum space, 

• setting e = in the convergent parts 
and performing a three-dimensional Fourier 
transform into coordinate space, 

• performing the po sums and in the end ana- 
lytically solving the remaining (hyperbolic) 
coordinate space integrals. 

The results for the purely bosonic sum-integrals 
are expressible in terms of a few mathematical 
constants, such as rationals. In 2, the Euler con- 
stant 7, and the derivatives of the Riemann zeta 
function C'(~l) ^-nd C'(~3). For the fermionic 
cases at finite chemical potential one in addition 
runs into the special functions C,'{x, y) = dx<^{x, y) 
and 4'i^) = r'(a;)/r(x) that occur in the combi- 
nations 

K(n,z) ^ C{-n,z) + {-ir^'C{-n,z*),{n) 
H(z) = *(z) + 1'(z*), (12) 



where z = \/2-i^i/{2iTT) and n e {0, 1, 2, 3}. For 
a somewhat detailed account of the asymptotic 
properties of these functions, see Ref. [12]. 

Proceeding further to four loops — one of the 
tasks required when approaching order g^ in the 
expansion of the pressure — the computations 
become considerably more involved. The sheer 
number of diagrams becomes so large that the 
automatization of the color and tensor algebras 
becomes imperative, and one in addition needs a 
more systematic and less time-consuming method 
for dealing with the remaining sum- integrals. 
This time the number of naive 'masters' will easily 
be 0(100), which one definitely needs to reduce 
in order for an analytic treatment to be possi- 
ble. The problem, however, is that most of the 
conventional methods for finding linear relations 
between integrals, most importantly the different 
integration- by-parts (IBP) routines, fail at finite 
temperature because of the discrete nature of the 
zeroth component of the integration momenta. 
New ideas are warmly welcomed. 

2.2. The soft scales 

The four-loop calculations one needs to per- 
form in the effective three-dimensional theories 
to obtain their contribution to the QCD pressure 
up to order g^lng — let alone g^ — are highly 
non-trivial already at the level of generating the 
diagrams and employing the IBP relations to find 
the master integrals (see Fig. 2). These computa- 
tions have, however, been described in some de- 
tail already in Ref. [13], so we will restrict the 
present treatment to cover merely the numerical 
evaluation of the four-loop masters performed in 
Ref. [14]. The strategy we chose to follow there 
was based on deriving linear difference equations 
in propagator powers for the integrals and solv- 
ing them numerically via factorial series, which 
was proposed and developed by S. Laporta in 
Ref. [15]. Below we comment on our implementa- 
tion of the algorithm only very briefly; for details, 
see Refs. [14,15]. 

Let [/ be a generic master integral. The basic 
idea behind the difference equation approach is 
to attach an arbitrary power x to one of its lines. 
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Figure 2. The one- (1), two- (0), three- (2) and 
four-loop (8) scalar master diagrams of EQCD 
[16]. The solid lines correspond to massive and 
dotted to massless scalar propagators. 



and by employing IBP identities in a systematic 
way to derive for it a linear equation of the form 



J2pji=^)u{x+j) 



F{x), 



(14) 



3=0 



U{x) 



1 



DID2...DN 



(13) 



where i? is a finite integer, the Pj's polynomials in 
X and d, and the right-hand side F{x) a function 
of simpler (already known) integrals. By substi- 
tuting a factorial series ansatz into Eq. (14) one 
can then derive recursion relations for the coef- 
ficients, and truncating the series at some finite 
Smax one finally obtains a result of desired accu- 
racy for the initial graph at some high x^ax ■ The 
last step is then simply to push the solution back 
to x = 1 by using the above difference equation, 
which however usually leads to a severe loss of ac- 
curacy in the results; hence a careful optimization 
of the ratio Smax/xmax is required. 

In our computations we employed the FORM 
program [17] to build the difference equations and 
recursion relations, and finally Mathematica to 
obtain the numerical solutions and expand them 
in powers of e. The first step utilized a slightly 
modified version of the IBP algorithm described 
in Ref. [13] to construct the equations and the 
routines introduced in Ref. [15] to solve them 
in terms of factorial series. The Mathematica 
part on the other hand consisted of numerically 
solving the recursion relations of the factorial se- 
ries coefficients, performing the actual summa- 
tion and finally implementing the push-back step. 
We used several tricks to increase the speed of 
the computations, some of which are explained 
in Ref. [14]. An essential feature in the project 
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Figure 3. The ring diagrams contributing to the 
QCD pressure at 0{g'^). 



was to attack the graphs in a specific order from 
the simplest one-loop case to the most compli- 
cated four-looper, which ensured that the inho- 
mogeneous terms of the difference equations were 
known at each step. 

3. SMALL T/n: EXPLICIT RESUMMA- 
TIONS 

To conclude, let us still briefly review the sta- 
tus of perturbation theory in the case of cold and 
dense quark matter, for which T/fi <C 1. It is ob- 
vious that in this limit the framework of dimen- 
sional reduction must cease to work, as the tem- 
perature no longer is the dominant energy scale 
in the problem. This can be explicitly verified 
in the exactly solvable large flavor number limit 
of the theory [18], where it is clearly seen that 
for values r//i < 1/10 the dimensionally reduced 
and exact results for the pressure start to deviate 
rapidly. In this region one is therefore forced to 
tackle the infrared problems — which again occur 
at three loops — through an explicit resummation 
of the IR divergent ring diagrams. These graphs 
are shown in Fig. 3, and one needs to add them to 
the sum of the known one-, two- and three-loop 
gluon-skeletons (2PI w.r.t. gluon lines) in order 
to render the whole result finite. 

The evaluation of the ring sum at arbitrary T 
and /i is a highly non-trivial problem even nu- 
merically. One simplification nevertheless occurs: 
since the vacuum part of the one-loop gluon po- 
larization tensor behaves as r^P^ in the IR limit. 



its contribution to the diagrams can be expanded 
as a power series in g^. This leads to a clean sepa- 
ration of all the singularities and renormalization 
scale dependence from the sum, and leaves an en- 
tirely finite quantity to be evaluated numerically. 
This task is being attacked at present [5]. 
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